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NON-INTERSECTING BROWNIAN BRIDGES AND THE LAGUERRE 

ORTHOGONAL ENSEMBLE 

GIA BAG NGUYEN AND DANIEL REMENIK 


Abstract. We show that the squared maximal height of the top path among N non¬ 
intersecting Brownian bridges starting and ending at the origin is distributed as the top 
eigenvalue of a random matrix drawn from the Laguerre Orthogonal Ensemble. This 
result can be thought of as a discrete version of K. Johansson’s result that the supremum 
of the Airy 2 process minus a parabola has the Tracy-Widom GOE distribution, and as 
such it provides an explanation for how this distribution arises in models belonging to 
the KPZ universality class with flat initial data. The result can be recast in terms of the 
probability that the top curve of the stationary Dyson Brownian motion hits an hyperbolic 
cosine barrier. 


1. Introduction and Main Results 

1.1. Motivation and background. The Kardar-Parisi-Zhang (KPZ) universality class 
describes a broad collection of models, including stochastic interface growth on a one¬ 
dimensional substrate, polymer chains directed in one dimension and fluctuating transver- 
sally in the other due to a random potential, driven lattice gas models, reaction-diffusion 
models in two-dimensional random media, and randomly forced Hamilton-Jacobi equa¬ 
tions. Although there is no precise definition of the KPZ universality class, it can be 
identified at the roughest level by its unusual scale of fluctuations (decorrelating at 
a spatial scale of The asymptotic distribution of the fluctuations, in the long time 

limit t oo, is conjectured to depend only on the initial (or boundary) condition imposed 
on each particular model. 

There are three special classes of initial data which stand out because of their scale 
invariance, usually referred to as curved^ flat and stationary. Based on exact computations 
for a few models which enjoy a special determinantal structure, the distribution of the 
asymptotic fluctuations in these three cases is known explicity. One of the most intriguing 
aspects of the KPZ universality class is that these limiting fluctuations are given in terms 
of objects coming from random matrix theory (RMT). This is particularly evident in the 
cases of curved and flat initial data: the asymptotic fluctuations are given, respectively, 
by the Tracy-Widom GUE and GOE distributions [TW94; TW96]. The first of these two 
distributions describes the asymptotic fluctuations of the largest eigenvalue of a random 
Hermitian matrix with Gaussian entries (the Gaussian Unitary Ensemble), while the sec¬ 
ond one is the analog in the real symmetric case (the Gaussian Orthogonal Ensemble); 
both will be introduced explicitly later on. For more background on this aspect of the 
KPZ universality class we refer the reader to the reviews [Corl2; QR14]; for some other 
perspectives we refer additionally to [Quail; BP14; QS15]. 

It is very natural in this context to wonder about what lies behind the connection between 
the KPZ class and RMT. Perhaps the most basic relationship one may seek is to find a 
model which lies in the KPZ universality class and which, at the same time, is naturally 
expressed as an object in RMT. As it turns out, in the case of the GUE (corresponding to 
curved initial data in the KPZ class) this can be achieved by considering a simple model: 
non-intersecting Brownian bridges (which we will introduce in detail in Section 1.2). This 
model is, on the one hand, one of the simplest and most studied models belonging to 
the KPZ class, while on the other hand it is equivalent to Dyson Brownian motion, a 
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process which describes the evolution of the eigenvalues of a GUE matrix whose entries 
undergo independent Ornstein-Uhlenbeck diffusions. A straightforward consequence of this 
equivalence is that the positions of the N non-intersecting Brownian paths at a single time 
are distributed as the eigenvalues of a GUE matrix of size N, and this leads directly to 
analog statements about their asymptotic fluctuations. Interestingly, the scope of this 
relationship extends also to looking at the entire paths of these processes. Eor instance, if 
one scales the top path of Dyson Brownian motion (or non-intersecting Brownian bridges) 
appropriately, then in the limit one obtains the Airy 2 process, which is known to describe 
the spatial fluctuations of models in the KPZ class with curved initial data. Beyond the 
basic relationship which we have just described, more recent developments in the area 
known as integrable probability have led to other, arguably deeper, ways of understanding 
the connection (see for instance [BP14; BG15]). 

The situation in the case of GOE, which corresponds to flat initial data in the KPZ class, 
is much less clear. In fact, essentially no results are available, and it has been a question 
of interest for several years now, both for probabilists and for physicists, to understand 
whether a relationship similar to the one available for the GUE case is available for GOE, 
or whether the appearance of the Tracy-Widom GOE distribution in the KPZ class is not 
much more than a coincidence. 

The fact that the GOE/flat link is much more difficult to understand is actually not 
surprising given that, as it is widely accepted, for most (if not all) models both in the 
KPZ class and in RMT, the GOE/flat case is considerably more difficult to analyze than 
the GUE/curved one. This is because many aspects of the integrability of these models 
which are present in the latter case, and lead to relatively simple exact formulas, are lost 
in the former one. It should be noted moreover that, in a certain sense, the GOE/flat 
connection is necessarily more tenuous than the GUE/curved one. In fact, if one considers 
the GOE version of Dyson Brownian motion then it is natural to expect (as conjectured 
in [BEPS07]), by analogy with the GUE case, that the top path would converge, under 
appropriate scaling, to the Airyi process, which is the analog of the Airy 2 process for 
models in the KPZ class with flat initial data. Nevertheless, [BPP08] provided convincing 
numerical evidence showing that this is not the case. 

The main goal of this article is to provide an explanation of how the GOE/flat link 
arises. We will achieve this by considering the model of non-intersecting Brownian bridges 
mentioned above but focusing now on a different quantity, namely the maximal height 
attained by the top path. Our main result will show that the distribution of the maximal 
height coincides with that of the largest singular value of a large rectangular matrix with 
Gaussian entries, or in other words, with the square root of the largest eigenvalue of a 
matrix from the Laguerre Orthogonal Ensemble, i.e. a real Wishart matrix. We remark 
that this identity will be established at the pre-asymptotic level (that is, for a finite number 
of paths and for a finite matrix), which is interesting in itself as we will explain in Section 
1.3. The connection with the GOE is established through the known RMT fact that the 
square root of the top eigenvalue of a real Wishart matrix converges under the right scaling 
to a Tracy-Widom GOE random variable. The way in which this result fits into the context 
of the KPZ universality class with flat initial data can be understood in terms of certain 
variational problems, and will be explained in Section 1.5. 

In the next two subsections we will change a bit our perspective to focus in more detail 
on the model of non-intersecting Brownian bridges, as well as on the Airy 2 process and on 
some previous results which relate it to the Tracy-Widom GOE distribution. 


1.2. Non-intersecting Brownian bridges. The model of non-intersecting Brownian 
bridges corresponds to considering a collection of N Brownian bridges {Bi{t),..., B]sf{t)), 
all starting from zero at time t = 0 and ending at zero at time t = 1, and conditioning 
them (in the sense of Doob) to not intersect in the region t G (0,1). We will always assume 
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that the paths are ordered so that Bi{t) < ■ ■ ■ < for t G (0,1). This model (which in 

the physics and combinatorics literatures is sometimes referred to as watermelons without 
a wall) and variants of it have been studied extensively in the last decade, see for instance 
[TW04; AM05; BS07; Fei08; KIK08; DKZll; FV12; Liel2; Johl3] among many others. The 
model can be thought of as a limit of non-intersecting random walks, which in the physics 
literature are known as vicious walkers, and were introduced by Fisher [Fis84] (under an 
additional conditioning on the walks staying positive) as a model for wetting and melting. 

The interest in studying systems of non-intersecting paths, both in the statistical physics 
and probability literatures, is due in large part to their intimate connection with RMT and 
the KPZ universality class. As an example, it has been shown that as the number of 
paths N —oo, and under proper scaling, several variants of these models converge to 
universal processes, such as the sine. Airy, Pearcey and tacnode processes. Universal here 
means that the same limiting processes arise for a wide class of other models (for more on 
this aspect see [Johl3; AFMIO] and references therein). The first two of these universal 
processes also arise naturally in RMT. For instance, and as we already mentioned, for fixed 
t G (0,1) the distribution of {Bi{t),... ,B]y{t)) coincides (modulo some scaling) with that 
of the eigenvalues of a random matrix drawn from the Gaussian Unitary Ensemble (GUE) 
and converge, under suitable scaling at the edge of the GUE spectrum, to the Airy point 
process. 

A particular aspect which has been subject of intense research has been the study of the 
maximal height attained by the highest path of a collection of non-intersecting paths. In 
the physics literature, [SMCRE08; Eei09; RSIO; RSll] obtain various expressions for the 
distribution of this maximum. As in the case of the limiting processes mentioned above, 
their main motivation lies in the computation of the asymptotic distribution in the N ^ oo 
limit, which for many different models is conjectured to be given by the Tracy-Widom GOE 
distribution. This was achieved in the physics literature using non-rigorous methods (see 
for instance [EMSll], which further establishes connections with Yang-Mills theory). Eor 
the case of non-intersecting Brownian motions on the half-line (with either absorbing or 
reflecting boundary condition at zero) this was rigorously proved by Liechty [Liel2]. 

In this paper we will focus on the distribution of the maximal height of a finite number 
of non-intersecting Brownian bridges. More precisely, for fixed N we are interested in the 
distribution of the random variable 

M-n = max BnU). (IT) 

te[o,i] 

As we already mentioned, under proper centering and scaling Mn should converge in 
distribution as ^ oo to a Tracy-Widom GOE random variable. The question in which 
we will be interested here is whether there is a finite N version of this result. Rather 
surprisingly, and as we mentioned already above, we will find that the answer is yes. But 
before stating the result, and in order to provide additional motivation (and in particular 
explain why this is in itself a natural question), let us discuss in some detail the GOE result 
in the N ^ oo regime. 

1.3. The Airy 2 process and GOE. The Airy 2 process A 2 was introduced by Prahofer 
and Spohn [PS02] in the study of the scaling limit of a discrete polynuclear growth (PNG) 
model. It is expected to govern the asymptotic spatial fluctuations in a wide variety of 
random growth models on a one-dimensional substrate with curved initial conditions, and 
the point-to-point free energies of directed random polymers in 1 -|- 1 dimensions. For its 
definition and a detailed discussion of its properties and relevance we refer the reader to 
[QR14]; let us just mention that the Airy 2 process is non-Markovian and stationary, with 
marginal distributions given by the Tracy-Widom GUE distribution. 

The Airy 2 process is also known to arise in the setting of (geometric) last passage per¬ 
colation. Here one considers a family {wii, j)}ij^z+ of independent geometric random 
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variables with parameter q (i.e. F{w{i,j) = k) = q{l — q)^ for k > 0) and let IItv be 
the collection of up-right paths of length N, that is, paths tt = (vro,..., 7 r„) such that 
TTj — TTi-i G {(1,0), (0,1)}. The point-to-point last passage time is defined, for M,N G 1,'^, 


by 

M-i-N 

= max w{'Ki), 

2 = 1 ) 

where the maximum is taken over all up-right paths connecting the origin to {M,N). 
Johansson [JohOO] proved that there are explicit constants ci and C 2 , depending only on q, 
such that P(LP““*( 1 V, A^) < ciiV-bC 2 iV^/V) —^ FovEir) as N ^ oo, with Fque the Tracy- 
Widom GUE distribution. Next one defines a process t e-?■ Hj\f{t) by linearly interpolating 
the values given by scaling LP°™*(A^, M) through the relation 


+ k,N - k) = ciN + C2N^/^Hnic3N-^/^k), (1.2) 


where C 3 is another explicit constant which depends only on q. Johansson [Joh03] went on 
to show that 

HNit)^A2{t)-t^ (1.3) 

in distribution, in the topology of uniform convergence on compact sets. On the other 
hand one can define the point-to-line last passage time by 

L'‘"®(iV)= max LP““VA^ + fe, iV - fe). (1.4) 

k=-N,...,N 

From the definition and Johansson’s result (1.3) it follows that 

C 2 -ciN]^ sup{^2(t) - 

teM 


in distribution. But it was known separately [BROl] that the quantity on the left converges 
in distribution to a Tracy-Widom GOE random variable, from which Johansson deduced 
in [Joh03] the remarkable fact that 

P Tmax {A2it) - t^) < = FGOE(4^'^^r), (1.5) 

Y j 

where Egoe denotes the Tracy-Widom GOE distribution (an explicit formula for Egoe 
will be given in Section 1.4). A more direct proof of (1.5) was given in [GQR13], based 
on formulas for the hitting probabilities for the Airy 2 process. This method has led to 
several other results about the Airy 2 and related processes (see e.g. [MFQR13] or the 
review [QR14]) and it is the one we will use in this paper in the context of non-intersecting 
Brownian bridges. 

The relation between the Airy 2 process and the study of AJat hes in the fact that, suitably 
rescaled, the top path of a collection of non-intersecting Brownian bridges converges to the 
Airy 2 process minus a parabola: 

2iVi/®(^R^(i(l+ lV-^/3^)) - Viv) -^A2{t)-t^ (1.6) 

in the sense of convergence in distribution in the topology of uniform convergence on 
compact sets. This result is well-known in the sense of convergence of finite-dimensional 
distributions; the stronger convergence stated here was proved in [GH14]. In view of this 
result, a similar argument as the one leading to (1.5) together with (1.5) itself gives the 
following: 


Theorem 1.1. 

lim pf2A^/®( max Bj^lt) — VTk) < r) = EGOE(4^^^r). (1.7) 

Af^oo V 4e[0,i] / 


It is this version of Johansson’s result (1.5) which provided the original motivation for our 
paper. We remark that, as a by-product of our results, we obtain a more direct derivation 
of (1.7). 
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1.4. GOE and LOE. In this section we will quickly introduce the two ensembles of ran¬ 
dom matrices which are most relevant to our results. The first one is the Gaussian Or¬ 
thogonal Ensemble (GOE). Let N'{a,b) denote a Gaussian random variable with mean a 
and variance b. An N x N GOE matrix is a symmetric matrix A such that Aij = AA(0,1) 
for i > j and An = AA(0,2), where all the Gaussian variables are independent (subject 
to the symmetry condition). The term orthogonal refers to the fact that the distribution 
of a GOE matrix is invariant under conjugation by an orthogonal matrix. The GOE can 
be regarded as the probability measure on the space of x real symmetric matrices A 
with density ) for some normalization constant Cn. The joint density of the 

eigenvalues (Ai,..., Aat) of a GOE matrix can be explicitly computed, and is given by 


1 

Zn 


N 

n 

2=1 


e 


l<i<j<N 


for some other normalization constant Z]\f. The weights appearing in this formula are 

the weights associated to the Hermite polynomials in the theory of orthogonal polynomials; 
for this reason, the Gaussian ensembles such as the GOE are sometimes also referred 
to as Hermite ensembles. The Wigner semicircle law [Wig55] states that the empirical 
eigenvalue density for the GOE has approximately a semicircle distribution on the interval 
[—2\/]V, 2y/N]. The fluctuations of the spectrum at its edges are of order and give 

rise to the Tracy-Widom GOE distribution: denoting by Agoe(A^) the largest eigenvalue 
of an At X At GOE matrix, we have [TW96] 

lim P(AGOE(At) < 2Vn + At-V6r-) = FcoEir) (1.8) 

N->-oo 

with 

-Pgoe(i’) = det(l - PoBj.Po)l 2 ('r), (1.9) 

where P^ denotes the projection onto the interval (r, oo) (i.e. Prf{x) = f{x)lx>r for 
/ G L2(M)), B, is the integral operator acting on L^(M) with kernel 


Br{x,y) = Ai{x+ y-Pr), 


(1.10) 


and Ai denotes the Airy function. The determinant in (1.9) means the Fredholm deter¬ 
minant on the Hilbert space L^(M). For the definition, properties and some background 
on Fredholm determinants, which can be thought of as the natural generalization of the 
usual determinant to infinite dimensional Hilbert spaces, we refer the reader to [QR14, 
Section 2]. We remark that (1.9) is not the original formula provided in [TW96] (which 
is written in terms of Painleve H transcendents instead of Fredholm determinants); this 
formula is essentially due to [Sas05], and was proved in [FS05]. Note also that one can 
choose a slightly different scaling (with the entries of a GOE matrix having variances N 
off the diagonal and 2N on the diagonal) so that the edge of the spectrum is at 2N and 
the fluctuations are of order which leads to a scaling in (1.8) of the same order as 

that in (1.2). 

We turn now to the Laguerre Orthogonal Ensemble (LOE). Let X be an n x A^ matrix 
whose entries are i.i.d. A/"(0,1), where we assume n > N. Then the matrix M = X'^X 
is said to be an X x X LOE matrix (often referred to also as a Wishart matrix). In 
applications to statistics, one thinks of the rows of X as containing n independent samples 
of an X-variate Gaussian population (with covariance matrix given by the identity), so that 
corresponds to the sample covariance matrix. The joint density of the eigenvalues of 
M is also explicit in this case, and is given by 

N 

n |A,-A,inAfe-^“/2, 

1<2<J<A^ 2=1 


1 

Zn 
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where the parameter a is defined to be a = {n—N—\)/2. The weights appearing in 

this case are the ones associated to the (generalized) Laguerre polynomials, which explains 
the name of this family of random matrices. By the Marcenko-Pastur law [MP67] the 
eigenvalues of M are concentrated on the interval [0,4A^]. Under our scaling, if a = 
{n — N — 1)/2 is fixed (and independent of N) then the fluctuations at the soft edge AN 
are of order and have the same limiting distribution as in the GOE case [JohOl]: 

denoting by Aloe(.A^) the largest eigenvalue of the LOE matrix, we have 

lim P(ALOE(iV) < 4iV + 2^/^N^I^r) = EGOE(r). (1.11) 

N^oo 

The scaling at the hard edge at the origin is different and gives rise to a different limit 
distribution, but we will not need it in this paper. 

In all that follows we will be interested exclusively in the case a = 0. 

1.5. Main results. We are ready now to state the main result of this paper. Let M be 
an LOE matrix with a = 0, that is, M = X'^X with X an (N + 1) x matrix with 
independent Af(0, 1) entries. Eor this choice of a we will denote by TLOE,Af the distribution 
function of the largest eigenvalue of M, 

pLOE,Af(i’) = IF’(Aloe(1V) < r). (1-12) 

Recall the definition in (1.1) of Mn as the maximum height of a collection of N non¬ 
intersecting Brownian bridges. 

Theorem 1.2. Let Bi[t) < ■■■ < B]\f{t) be a collection of non-intersecting Brownian 
bridges as above. Then for all r > 0 we have 

p( max y/2BN{t) < r ] = FLOE,A(2r^). (1-13) 

Vie [0,1] ) 

In other words, 4A4^ is distributed as the largest eigenvalue of an LOE matrix or, alter¬ 
natively, 2Xijsi is distributed as the largest singular value of the {N + 1) x matrix X 
introduced above. 

Let us quickly verify that the scaling in this result is consistent with the one in Theorem 
1.1 and (1.11). Theorem 1.2 says that ^ Aloe(A7)/4. By (1.11), this implies that 

Mn = ^Ai + 2-2/3 aU/3^goe + o(AiV3) = ViV + 2-^I^N-^I^C,go^ + o{N-^l^), 

where Cgoe is a Tracy-Widom GOE random variable. In other words, 2Ai^/®(A4Ar~V^) = 
4 “^/^(^goe + o(l)) which is exactly the content of Theorem 1.1. In particular, Theorem 1.1 
follows as a corollary of (1.11) and (1.13). 

We take a brief pause now and go back to an issue left open at the end of Section 1.1, 
which is the question of why Mat should be interpreted as a flat initial data object in the 
KPZ universality class. In a way, the convergence of Mw a Tracy-Widom GOE random 
variable should be taken, in itself, as enough evidence of this fact. But the connection 
goes a bit further, and can be understood in terms of certain variational formulas. For 
example, in the context of last passage percolation (LPP), the point-to-line last passage 
times (1.4) leading to GOE fluctuations are defined in terms of the maximum of point- 
to-point last passage times (1.2), which in turn lead to GUE fluctuations. The parallel 
with (1.13) is direct. The exact same relationship can be established at the level of many 
other polymer models (of which LPP is a zero-temperature version), and at the level of 
the totally asymmetric exclusion process (which can be mapped to LPP). 

This straightforward way of expressing flat initial data quantities in terms of their curved 
initial data analogues is not as explicit in the case of some other models, such as the partially 
asymmetric exclusion process, which have less (or at least a more complicated) algebraic 
structure, but it is interesting to note that it does hold at the level of another of the 
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most important members of the KPZ universality class, the KPZ equation. Without going 
into much detail, the KPZ equation can be understood by studying the stochastic heat 
equation (SHE), which is linear. It turns out that the flat initial data for the KPZ equation 
corresponds to starting the SHE with initial condition Zq = 1, and thus by linearity the 
flat solution can be obtained by convolving the constant function 1 with the solution of 
the SHE starting with Zq = 6o, which corresponds to curved initial data. Note that the 
relationship in this last case is not written directly in terms of a variational problem as 
described before, but one can check that (at least conjecturally, by essentially appealing 
to a version of Laplace’s method) one recovers a variational problem as time t —)• oo. Eor 
much more on this see [QR14; CQR15; QR15]. 

Coming back to the description of our main results. Theorem 1.2 is equivalent to a state¬ 
ment about the probability that the top path of Dyson Brownian motion hits an hyperbolic 
cosine barrier, and it is that version of the result which we will prove. Consider an N x N 
random matrix drawn from the Gaussian Unitary Ensemble, that is, a (complex-valued) 
Hermitian matrix A such that Aij = A/"(0,1/4) -|- iA/"(0,1/4) for i > j and An = A/"(0,1/2), 
where all the Gaussian variables are independent (subject to the Hermitian condition). 
Note that (for later convenience) we have chosen a somewhat different scaling for the 
Gaussian variables here compared with our definition of the COE matrices. Now suppose 
that we let A evolve by letting each Gaussian variable in the construction diffuse according 
to independent copies of the Ornstein-Uhlenbeck process Xt defined as the solution of the 
SDE 


dXt = -Xtdt + a dWt, 


where Wt is a standard Brownian motion and ^ off-diagonal entries and cr = 1 

on the diagonal. We write the eigenvalues of this matrix at time t as (Ai(t),... 
with Xi{t) increasing with i. This eigenvalue diffusion is known as the stationary (GUE) 
Dyson Brownian motion, and it defines an ensemble of almost surely non-intersecting curves 
indexed by M. It can alternatively be written as the solution of a certain A^-dimensional 
SDE, and it is not hard to check that it is stationary, with marginals at any time t given 
by the eigenvalue distribution of an x GUE matrix. 


Theorem 1.3. Let (Xi{t),..., Xisf{t)) be the stationary Dyson Brownian motion defined 
above and let FLOE,Af be defined as in (1.12). Then 


F{XN{t) < rcosh(t) Vt G M) = ELOE,Af(2r^). (1-14) 

The equivalence between the two results is due to the fact that non-intersecting Brownian 
bridges can be mapped into the stationary Dyson Brownian motion in such a way that the 
probabilities on the left-hand side of (1-13) and (1.14) coincide. We will explain this in 
more detail in Section 2. 

The proof of Theorem 1.3 has two steps. The first one consists in obtaining an ex¬ 
plicit formula for the probability on the left-hand side of (1.14). By the mapping between 
non-intersecting Brownian bridges and the stationary Dyson Brownian motion alluded to 
above, this is equivalent to finding a formula for the distribution of Mn- As we already 
mentioned, there are formulas in the literature for the distribution of the maximal height 
of several models related to non-intersecting Brownian bridges, which can be obtained 
through a direct application of the Karlin-McGregor/Lindstrom-Gessel-Viennot formula 
[KM59; Lin73; GV85]. Eor completeness, let us state the formula in the case of AIat (see 
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[SMCRF08])!; 


^2N 


P(A^7V <r) = 




dye 




det 


Vi ^cos(yi + ^) 


N 


J i,j=l 


(1.15) 

By using the Cauchy-Binet identity, the right-hand side can be turned into a single N x N 
determinant with entries involving Hermite polynomials, see (102)“(103) in [RSll]. The 
resulting formula is reminiscent of some of the formulas we will obtain below, see (3.6) 
together with (1.17), but it is not clear how to use it directly to obtain a proof of Theorem 
1.1 (nor of (1.13)). Moreover, as we will explain next, while the structure of the Fredholm 
determinant formula for the distribution of which we will obtain in this paper (see 
Proposition 1.4) makes very apparent a connection with Johansson’s result (1.5) — this 
was an important clue for us in the discovery of (1-13) — from the formula appearing in 
[RSll] such a connection is not at all clear. It is worth mentioning that in the case of 
Brownian excursions, for which the analog of (1.15) turns out to be slightly simpler, the 
analog of Theorem 1.1 (with the same limit) was proved by Liechty [Liel2] by appealing 
to a Riemann-Hilbert analysis of a certain system of discrete orthogonal polynomials. 


Here we follow a different strategy, leading to an arguably simpler formula which also 
has some intrinsic interest. Working at the level of Dyson Brownian motion, we appeal to 
a result of [BCR15] in order to obtain an expression for < rcosh(t) Vt € [—L,L]), 

for fixed L > 0 in terms of the Fredholm determinant of what they call a “path-integral 
kernel”. This path-integral kernel can be expressed in terms of the solution to a boundary 
value PDF, which we then solve explicitly. Taking L —)■ oo in the resulting formula leads 
to the following result. Let (pn be the harmonic oscillator functions (which we will refer 
to as Hermite functions), defined by (pn{x) = e~^ ^‘^pn{x), with pn the n-th normalized 
Hermite polynomial (i.e., so that ||¥?||2 = 1), and define the Hermite kernel as 


N-l 

Knerm, N{x,y) = ^ pn{x)pn{y)- (1T6) 

n=0 

We introduce also the reflection operator Qj. on L^(M), defined by 


Qrf{x) = f{2r - x). 


Proposition 1.4. For any r > 0, 

P(AAr(t) < rcosh(t) Vt € M) = det(l - KHerm,N£>rKHerm,Af)A 2 (M) • (1-17) 

The same formula holds for P(maxjg[o,i] V^B^it) < r). 

This result will be proved in Section 2. 

The expression on the right-hand side of (1.17) is a close analog of the formula for Fgoe 
appearing in (1.9). To see this we introduce the Airy kernel, defined as 

poo 

KAi{x,y) = / dXAi{x + X)Ai{y + X). 

Jo 

This kernel is closely related to GUE, as it is the limiting correlation kernel of the GUE 
eigenvalues near the edge of the spectrum. It is related to the Tracy-Widom GOE distri¬ 
bution because of the identity dXAi{a + A) Ai(6 — A) = 2“^/^ Ai(2“^/^(o J- b)), which 
(since Kai = BoPqBo, with Bq defined in (1.10)) implies that 

KAi^'rKAi = BoPoB^PoBo (1-18) 

^This formula was derived in [SMCRF08] using path-integral techniques. Although we are not aware 
of a derivation in the literature based on the Karlin-McGregor formula, for the case of non-intersecting 
Brownian excursions (corresponding to imposing an absorbing boundary at zero) the analog formula, also 
derived in [SMCRF08], was rederived in this way in [KT07]. 






NON-INTERSECTING BROWNIAN BRIDGES AND THE LAGUERRE ORTHOGONAL ENSEMBLE 9 


with Br(x,y) = 2“^/^ Ai(2“^/^(x + y + 2r)). Since Bq = I (this identity is related to the fact 
that the family of functions {Ai(x + constitutes a generalized eigenbasis of L^(M)), 

the cyclic property of the determinant and (1.9) allow us to conclude that 

.fGOE(4^/^r) = det(l — KAi£»rKAi)L 2 (K) . (1-19) 

We point out that there does not appear to be a direct analog of (1.18) for Knerm.Af 
(although one can obtain explicit formulas for KHerm, 7 V 0 rKHerm, 7 V involving no integrals, 
see for instance (A.4) and (A.6)). 

We can actually push the analogy between (1.17) and (1.19) a bit further and use it to 
provide a simple proof of Theorem 1.1. Indeed, a simple scaling argument on the right- 

hand side of (1.17) leads to P^2A7^/®(A4Ar — '/N) < rj = det^l — KHerm,Ar^rKHerm,N^ with 

KHerm,Af(3:, U) = Knerm.Af(«^AfiC + , KnU + V^), where Kjsf = On the 

other hand, it is well known that KHerm,Ar converges to Kai as A" —>■ oo, where the conver¬ 
gence is strong enough to imply the convergence of the associated Fredholm determinants. 
In view of (1.19), and omitting the details, this implies Theorem I.l. 

A related observation is that, in a sense, Proposition 1.4 serves as a generalization of 
Johansson’s result for the Airy 2 process, (1.5). In fact, the scaling argument used in the 
last paragraph leads to det(l — KHerm,Af^’rKHerm,Ar ) = P(AAr(i) < (/tAfi’T V2A) cosh(t) Vt G 
M). On the other hand, it is known that AAr(i) = ~ V^N) converges 

to A 2 {t) (this is just a restatement of (1.6) in view of the mapping between the two 
models). If we knew that the convergence is strong enough to imply the convergence of 

P^AAr(t) < a Vt G with some control on a, then (1.5) would follow, because by the 

argument sketched in the last paragraph the determinant would go to FGOE(4^'^^r), while 
[(kati’ + V2A) cosh(A“^/^t) — -\/2iVj = r -|- -|- 

The second step in the proof of Theorem 1.3 consists in showing that the right-hand 
side of (1.17), i.e. det(l - KHerm,Af£'rKHerm,N)A2(R), equals FLOE,iv(2r^). This is proved in 
Section 3. We remark that, together with the preceding discussion, this identity provides 
an alternative proof of the result of [JohOl] in the case a = 0. 


2. Hitting probabilities for Dyson Brownian motion 

Recall the stationary Dyson Brownian motion introduced in Section 1.5. As we men¬ 
tioned, this model is intimately related to non-intersecting Brownian bridges. The basic 
relation is that if one considers the non-stationary version of Dyson Brownian motion 
(where the Gaussian variables making up the entries of a GUE matrix evolve according to 
a plain Brownian motion), then the dynamics of the eigenvalues of this evolving matrix 
coincide with those of a collection of Brownian motions conditioned to never intersect. 
The analogous relation in our setting goes through a time-change, and is given explicitly 
in [TW07, Section 2.2.2]: if Bi{t) < ••• < t G [0,1], are non-intersecting Brown¬ 

ian bridges and Ai(t) < • • • < AAr(i), t G M, are dehned as a stationary Dyson Brownian 
motion, then 


Ai(ilog(t/(l -t)))),^^ ^ 

as processes defined for t G [0,1]. Ghanging variables 1 1 —> e^^/(l -|- e^®) leads to 


max B]\f{t) = 
tG[0,l] 


max 

tG[0,l] 


^2t{l - t) Xn{^ log{t /(1 - t))) =sup 

seR v2cosh(s 
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which shows that Theorems 1.2 and 1.3 are equivalent^. The rest of this section will thus 
be devoted to computing P(Aw(t) < rcosh(t) Vt G M). 


2.1. Path-integral kernel. The finite-dimensional distributions of the stationary (GUE) 
Dyson Brownian motion are classically expressed through a Fredholm determinant in terms 
of the extended Hermite kernel ^ 


Kext j T.nJ ^^^n{x)^n{y) if S > t, 

ifxt 

where (pn{x) = e~^^^^pn{x) and pn is the n-th normalized Hermite polynomial. Explicitly, 
if —oo < ti < t 2 < ■ ■ ■ < tn < oo and n,..., G M, then 


P(Aiv(ti) < Vj, j = 1, . . . ,n) = det(l - , (2.1) 

where we have counting measure on {ti,... ,tn} and Lebesgue measure on M, and f is 
dehned on {ti,... , tn} x M by 

i[tj,x) '^xe{rj,oo) 


(for more details see [TW07]). 

The first step in our derivation is to obtain a formula for the probability that A 7 v(t) 
stays below rcosh(t) on a hnite interval [—L,L]. To that end, we need to consider a finite 
mesh ti < • • • < of [—L, L], let r* = r cosh(tj), and then take a limit of the corresponding 
probabilities as given in (2.1) as the mesh size goes to zero. But these probabilities become 
increasingly cumbersome as n increases, due to the n-dependence in the L? space on which 
the operators act. The way to overcome this problem is to first manipulate the right-hand 
side of (2.1) into a Fredholm determinant of some other kernel acting on L^(]R). Such a 
formula was first stated, in the context of the Airy 2 process, in [PS02] (see also [PSll]), 
and the resulting formula was used in [CQR13] to obtain a formula for the probability 
that A 2 {t) stays below a given function g[t) on a finite interval. Later on, the procedure 
that converts the extended kernel formula into a formula with a Fredholm determinant 
acting on L^(M) was generalized in [BCR15] (see also [QR13]) to a wide class of processes 
that includes the stationary Dyson Brownian motion, and from the resulting formula they 
obtained a continuum statistics formula for Dyson Brownian motion in a similar way as in 
[CQR13]. In order to state the formula we need to introduce some operators. 

First, recall the definition of the Hermite kernel KHerm.iV) given in (1.16), and note that 
^iferm y) ~ KHerm,Af(a^j 2/) for any t. Next we introduce the differential operator 

D = -i(A-x2 + l) 

(A is the Laplacian on M). D and KHerm.w are related: Dipn = npn, so that KHerm.w 
is the projection operator onto the space span{(/?o,..., associated to the hrst N 

eigenvalues of D. In particular, even though is well-defined in general only for t < 0, 
e*^KHerm,Af IS Well defined for all t, and its integral kernel is given by 

N-l 

e*°KHerm,w(a:,2/) = ^ e^'^ipn{x)ipn{y)- (2.2) 

n=0 

Now hx ii < £2 and consider a function g G H^{[ii,£ 2 ]) (he. both g and its derivative 
are in L^(M)). We introduce an operator 0^^^ acting on L^(M) as follows: 0^^^ ^^j/(x) = 


similar argument, together with the fact [TW04] that — y/2N) converges to 

A 2 {t) in the sense of finite-dimensional distributions, provides a justification for a version of (1.6) in this 
weaker sense. 
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u{(. 2 ,x), where ix(£ 2 , •) is the solution at time £2 of the boundary value problem 

dtu + Du = 0 for a: < g{t), t G (t'l,£ 2 ) 
u{£i,x) = (2.3) 

u{t,x) = 0 for a: > g{t). 

Proposition 2.1 ([BCR15]). For any £i < £.2 and g G H^{[£\,£ 2 ]) we have 

P(AAr(t) < g{t) Vt G [ 4 ,^ 2 ]) = det(^l - Knerm.Af + • (2.4) 


See [BCR15, Proposition 4.3 and Remark 4.4] for more details. Here, and in the rest of 
this section, the Fredholm determinant is computed on the Hilbert space L^(R). 

In order to make use of (2.4) we need a formula for 0^^^ By the linearity of (2.3), 
®f^i h] integral operator with kernel given by solving the boundary value problem 

with / replaced by a delta function. The next result gives a probabilistic representation 
for the integral kernel of 0^^^ . 

Proposition 2.2. Let a = /3 = and denote by ^^^{x,y) the integral kernel 


ely^^{x,y) = e-^(y 




i/47r(/3 - a) 


b{a)=e^'i^x, S(/3)=e^2y 


< \/4tc/(4log(4t)) Vt G [a,/3]^ , (2.5) 


where the probability is computed with respect to a Brownian bridge b{t) from e^^x at time 
a to e^^y at time (3 and with diffusion coefficient 2. 


Proof. Let u{t, x) be the solution to the boundary value PDF (2.3) and consider the trans¬ 
formation u{t,x) = e* z) with r = je^^, z = e^x. It is not hard to check then that 

v{t, z) satisfies the following boundary value problem associated to the heat equation: 

drV — d‘^v = 0 for z <-\/4r (7(log(4r)/2), T£{a,f3) 

v{a,z) = 

v{t,z) = 0 for 2 > \/4rg'(log(4r)/2), 

where a = /3 = This boundary value PDF can be solved explicitly in terms 

of Brownian motion by using the Feynman-Kac formula: letting b{s) denote a Brownian 
bridge with diffusion coefficient 2, we have 


yVIag(log(4a}/2) 2 //o ^ n/o / X \ 

v{/3, z)= dx e-V(8a)-log(4a)/2^ / ^ 

7-00 ^/47r{f3 - a) 

■^h{a)=xMfi)=z ^ V^5(log(4r)/2) on [a,/3]) . 

Now using the fact that u{£ 2 ,y) = e^^/4, e^^y) and recalling that a = ^e^^^ we 
immediately obtain 

fePgiP) ^ 2 ^ g-(x-e^2j/)2/(4(/3-a)) 

u{£2,y)= dxe^y 2 ® *+2 1 - — f{e^x) 

J -00 ^f4:7^{/3 — a) 

■^b{a)=x,m=P 2 y {Hr) < V^g{log{4T)/2) on [a,/3]) . 

Changing variables x i-a in the integral, the formula for 0^^^ (x, y) readily follows. □ 
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2.2. Hyperbolic cosine barrier. Observe now the key fact that, in our case g{t) = 
rcosh(t), the probability appearing in (2.5) is reduced to the probability of a Brownian 
bridge staying below the linear function 2rt + ^r, which can be computed explicitly. In fact, 
assuming that x < e~^^{2ra + r/2) = rcosh(^i) and y < e~^‘^{2rj3 + r/2) = rcosh(£ 2 ) (note 
that the probability below is obviously zero if either condition is not met), the Cameron- 
Martin-Girsanov formula yields 


b(a)=e^l X 
b(/3)=e^2y 

= 1 — e 


< 2rt + on [a, 


—r(e^2 y—e^l x)+r'^ (fi—a) 


(e^2 y —2r/3 —e^l x+2ra)‘^ 
g 4(/3-c«) 


(e‘2y-Jlx)2 b{a)=e‘lx-2ra \ 

e 4(13-0,) b{p)=e‘2y-2rp ^ ^ 


max b{t) > hr ] . 


The last probability can be computed easily using the reflection principle, and it equals 
g-(e^ia;- 2 ra-r/ 2 )(e^ 2 y- 2 r/ 3 -r/ 2 )/(/ 3 -o) instance page 67 in [BS02]). Putting this back 

in our formula (2.5) for with g{t) = rcosh(t), which for simplicity we will denote 

(r) 

[b,^ 2 ]’ 


from now on as « i, gives 


^x<rcosh{ii),y<rcosh{£ 2 ) ^ 


1 


[b,f 2 ] 

X (e 


i/47r(/3 - a) 


— (e^lx—e^ 2 y) 2 /( 4 (, 3 _Q,)) r{e^2y—e^3x)+r'^{fi—a)—{e^3x+e^2y—2r{a+j3)—r)^/{A{l3—a)) 


( 2 . 6 ) 


The above expression splits into two terms. Note that if we disregard the indicator function, 
then by the above derivation the first term corresponds to the solution of (2.3) with g = oo, 
and thus it is nothing but {x,y). As a consequence, we deduce that 


0h) _ p 

[Hi,£ 2 ] '^rcos 


(p~{^2—£\)^ _ p(^) p 

ih(£i) (^6 ^[^ 1 / 2 ]/ ^»'cosh(£2)’ 


(2.7) 


(r) 

where Paf{x) = (I — Pa)f{x) = f{x)lx<a and RL% 1 is the reflection term 


[^l:-^ 2 ] 


R 


li.’fcA'!') = 71 ^'’ 


e^ix)+r^(/3—a)—(e^ix+e^2y—2r(a+/9)—r)^/(4(/3—a)) 


( 2 . 8 ) 


and, we recall, a = /3 = 

Now we set —£i = £2 = L, so that by Proposition 2.1 we have 

F{XN{t) < rcosh(t) Vt G M) = ^i^det^^l - KHerm,Ar + • 

Using now the cyclic property of the Fredholm determinant and the identities e^'^*^KHerm,Af = 
(e RHerm,Af) and e KHerm,AfC RHerm,A' ~ ^ RHerm,Af6 RHermjAl ~ RHerm,Af (which 
follow directly from (2.2) and the orthonormality of the (pns) we may rewrite the last iden¬ 
tity as 

IF’(A7v(i) < r-cosh(t) Vt G M) = ^li^det(^l - KHerm,A + e^^Kuerm,if ■ 

(2.9) 

Note that s 1 —> is a semigroup, so that o]®[o L]’ ™ view 

of (2.7) we may write 


e|:7i, = P„„h(L)(e-“ - - R}„T)P-o. h(L) 


-LD dT) 


-LD d('-) 


fr’) 

Following [CQR13], we decompose 0[_j^ in the following way 


iP plr-J -N qI*- 


?h) 


-LD oh) 


ih) 


( 2 . 10 ) 
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where 




- R 


(r) 

[-L,0] 




Prcosh(L) 


-LD 



f^[ 0 ,L]) ^rcosh(L)- (2-11) 


(r) 

The idea is that is an error term which goes to 0 as L —?■ oo. This is the content of 
the next result, whose proof we defer to Appendix B: 


Lemma 2.3. Assume r > 0. 
norm. 


Then := e^'^KHerm,lvf^£^e'^^KHerm,Ar -1 0 in trace 

’ ’ L^OO 


Since the mapping A i—)■ det(l + A) is continuous with respect to the trace norm, the 
lemma together with (2.9) and (2.10) show that if 


A : = 


lim 

L —^CXD 


KHerm,Af “ Knerm.Af (« 





( 2 . 12 ) 


exists in the trace class topology, then 

P {Xisf{t) < rcosh(t) Vt € M) = det(l — A). (2-13) 

But, as we will see next, the operator inside the brackets in (2.12) in fact does not depend 
on L (the analogous property was proved in [CQR13] in the setting of the Airy 2 process). 
The key step is the following result: 


Lemma 2.4. For all L > 0, 

i/D (^) (^) TjD 

C I^Herm,7V q] ~ ^^Herm,Af Qr and ^[0,L]® ^Herm,7V ~ ^^Herm,N j 

where Qr is the reflection operator Qrf{x) = f{2r — x). 


Using this lemma and the fact that e'^’^Knerm.Af e = KHerm.iv we get 

Knerm.Af “ ^Rerrn.N KHerm,Ar 

— RHerrOjAf RHerm,A'(l i?r)Pr'(l i?r) RHerm,Af — RHerm,7V ^’r RHerm, A"; 

where the second equality follows from the identity (I — ^?r)Pr(l ~ ^r) = I ~ i?r- In view of 
(2.12) and (2.13), this yields Proposition 1.4. All that is left to prove then is Lemma 2.4. 


Proof of Lemma 2.4- We will only provide the proof of the first formula, the second one is 
very similar. Using (2.8) we write the kernel of the operator R^^. o] 




V^7r(l-e 2r) 


^—ax^+byX+Cy 


with 


a = - 


, _ 2e ^{2r-y) 

~ l-e-2^ 


and 


2(l-e-2r) 


2(l-e-2r) ! 

This formula together with (2.2) and the contour integral representation of the Hermite 
function (^n(a^), 

r,2tX — t^ 


(/7„(x) = (2”n!v^) ^ d,t- 


tn+l 


(where the contour of integration encircles the origin), gives us 




N-1 


(r) 


n=0 


^yTT(l-e 


N-1 I „ -P r 
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The 2 : integral is just a Gaussian integral, and computing it the last expression becomes 


N-l , 

e^>4x)(2"n!0F)-'/2^ f dt- 

n=0 


^-e-^^t^+2e-^t[2r-y)-{2r-yY/2 


^n+l 


N-l , 

j) dt- 

n=0 


j-t2-H2t(2r-j/)-(2r-j/)2/2 
pi+l^Ln ' 


where we have performed the change of variables t te^. The last integral and its prefac¬ 
tor are nothing but ipn{2r - y), so this yields e-^°KHerm,AfR[l.i^o](^’ = KHerm,Af(a;, 2r - y) 
as needed. □ 


3. Connection with LOE 


This section is devoted to the proof of the following result: 
Proposition 3.1. For r > 0, 

det(l — KHerm,AferKHerm,Af)i 2 (R) = ELOE,Af (2r^). 


Together with Proposition 1.4, this proposition implies Theorem 1.3. 

Let us start by introducing an explicit formula for FloEjN- To that end, we will utilize 
the ensemble A(l) < A(2) < • • • < A(iV) obtained as the result of superimposing the 
eigenvalues of two independent copies of our LOE matrices, writing them in increasing 
order, and then keeping only the even labelled coordinates (i.e. keeping the largest, 3’^'^ 
largest, 5**^ largest, and so on). Observe that if Aloe(-A^) denotes the largest eigenvalue of 
an LOE matrix as in Section 1.4, then 

P(ALOE(iV) < 2r2)2 = P(A(iV) < 2r2). (3.1) 

The advantage of this representation is that the superimposed ensemble is a 

determinantal process with a simple correlation kernel KLague, 7 V (see [FR04]). The kernel 
KLague,Af IS given as follows. For n € N, introduce the Laguerre function 

ifnix) = e"*/^L„(x), (3.2) 

where Ln{x) is the n-th normalized Laguerre polynomial (so that HV’nlb = 1); and then 
dehne the Laguerre kernel as 


N-l 

KLague,Ar(T,I/) = ^ 'llln{x)lpn{y) ■ 
n=0 


Then 


d fy 

RLague,A'(®) y) ~ J KLague,Af(3:, u). 


The determinantal structure of the superimposed ensemble leads directly to a formula for 
the distribution of X{N) (see [Joh03] or [QR14, (1.36)]): 


P(A(A^) < 2 r^) = detfl - P2r2KLague,AfP2r2) • ( 3 - 3 ) 


Observe that KLague.N is a finite rank operator, and thus the last determinant can be 
represented as the determinant of a finite matrix. More precisely, if we factor our op¬ 
erator as Knerm.Af = K 1 K 2 with Ki : £^({0,... ,N — 1}) —L^(]R) and K 2 : L^(M) —> 
.£^({0,... , — 1}) defined by the kernels Ki(x,n) = —ipn{x) and K 2 (n, y) = Jq duifniu), 
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then the cyclic property of the Fredholm determinant implies that the determinant in (3.3) 
equals det(l — K 2 K 1 ), so that 


^ r /*oo px 

det(l - P 2 r- 2 KLague,AfP 2 r 2 )i 2 (R) = det 6jk + J ^dxip'jix) J dui/jhiu) 


N-1 
J j,k=0 


= det 


1*00 

6jk- / dx'il;j{x)'ipk{x) - 'ipj{2r'^) / duipki 
J 2 r 2 Jo 


-1 N-1 

\ j,fc=o’ 

where in the second equality we have integrated by parts. Defining now a symmetric 


matrix^ L G 


Ljk — 


tiNxN 


and two column vectors Ri,R 2 G by 


> p2r^ 

dx'ipj{x)'ijjk{x), {Ri)j = ipj{2r‘^) and (i? 2 )j = / du'ijjj{u), (3.4) 

Jo 


/2r2 Jo 

for j,k G {0,... , N — 1}, we deduce by the last identity, (3.1) and (3.3), that 

■PLOE,Ar(2r^)^ = det(/ — L — Ri® R 2 ). 


(3.5) 


Similarly, we have a version of (1.17) in terms of the determinant of a finite matrix (which 
can be obtained by conjugating the kernel inside the Fredholm determinant in (1.17) by 
the operator G: L^(M) —> ^^({0 ,... ,N — 1}) with kernel G{n,x) = (fn(x)): 

det(l KHerm,Af l^Herm,Af) — det(/ 77), (3'6) 

where the symmetric matrix ff has entries given by 

Hjk= / dxipj{x)tpk{2r - x). (3.7) 

JR 

Therefore, and in view of (3.5), we see that, in order to prove Proposition 3.1, we have to 
establish that 

det{I - Hf = det{I - L-Ri® R 2 ). (3.8) 


At this point the main difficulty in proving (3.8) lies in the fact that the two sides of 
the identity are given in terms of objects related to two different families of orthogonal 
polynomials, which makes it hard to relate one to the other. So the first step in our proof 
of the identity consists in replacing the matrix H on the left-hand side by a matrix defined 
in terms of Laguerre polynomials. 

To this end, let us introduce the following N x N (real) matrix H: 

Hij = (-1)^ (V’*-hi-Af(2r^) - 'ilJi+j-N+ii‘2r‘^)) for i,j G {0,... ,iV - 1}. (3.9) 

Here ifin is the Laguerre function introduced in (3.2) for n > 0, while we set = 0 
for n < 0. Note in particular that H is zero above the anti-diagonal (i.e. Hij = 0 if 
i + j < N — 1). This matrix will play a key role in the proof. As we will see in the next 
lemma, H is conjugate to 77, so that det(7 — 77) = det(7 — 77). Moreover, we will see that 
the matrices L and 7?i ® R 2 are also intimately related to 77. In order to state the lemma 
we introduce a matrix Q G and two column vectors u,v G by 

u= (- 1 )^- 4 , Vi = {-iy2 for i = 0,...,A^-l, 


Qij 


0 for i < j, 
< —2r for i = j, 
—4r for i > j, 


i,j = 0,... ,N -1 


(3.10) 


(here 1 denotes the constant vector with 1 in each entry). 

Note that the matrices and vectors introduced in this section are always indexed by 
{0,... , TV — 1}, and they generally depend on N and r; we have omitted this dependence 
from the notation for simplicity. 


^As a notational guide, note that while we have used sans-serif fonts to denote operators acting on a 
Hilbert space and their associated kernels, we are using regular fonts to denote (finite) matrices. 
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Lemma 3.2. Let H, H, L, Ri, R 2 , Q, u and v be defined as in (3.4), (3.7), (3.9) and 
(3.10). Then the following properties hold: 

(i) H is conjugate to H, i.e. there exists an invertible matrix S G such that 

H = SHS-fi 

(ii) H‘^ = L- 

(iii) Ri = Hu and R 2 = (I — H)v. 

(iv) ^H = QH. ^ ^ 

(v) |:(/ + H)-^ = {I- + (/ - + H)-^, where E = ArHu ® u. 

This lemma contains all the key identities which will be needed in the proof of (3.8). Let 
us thus postpone the proof of the lemma until the end of this section and proceed directly 
to the proof of the main result of this section. 


Proof of Proposition 3.1. As we already explained, all we need to do is prove (3.8). The 
structure of the proof is inspired in that of the proof of (1.9) in [FS05]. Note that both 
sides of (3.8) are zero if r = 0 (this is equivalent to the fact that both sides of (3.1) vanish 
when r = 0, which is clear). Therefore we will assume throughout this proof that r > 0, 
which for similar reasons implies that both det(/ —77) and det(7 —L —7ii(8)7?2) are strictly 
positive. 

We start by using (ii) and (iii) in Lemma 3.2 to rewrite the determinant on the left-hand 
side of (3.8) as 

det(7 - L - 7?i ® 7 ^ 2 ) = det(7 - - Huv'^{I - H)) 

= det(7 + H) det (7 - (7 + Uy^Huv'^) det(7 - H) 

= det(7 - H) det(7 + H){l- {u, (7 + H)-^Hv )), 

where in the third equality we used the fact that (7 -|- H)~^Huv^ is rank one. By Lemma 
3.2(i), we have det(7 — 77) = det(7 — 77), and thus (3.8) will follow if we prove that 

det(7-.H') = dei{I + H){l - {ufil + HY^Hv)). (3.11) 

Note that, by the discussion in the last paragraph, since r > 0, the left-hand side and the 
two factors on the right-hand side are strictly positive. 

Consider the second factor on the right-hand side of (3.11). Since {u, u) = 0 if is even 
and (u, u) = 2 if is odd, we can write 1 = (u, v) + (—l)'^, so that 

1 - (u, (7 + H)-^Hv) = (-1)^ + (u, v) - {u, (7 + H)-^Hv) = (-1)^ + {u, (7 + HyY). 


Taking now logarithm on both sides we see that (3.11) is equivalent to 

log det(7 — 77) = log det(7 -|- 77) -|- log((—1)^ -|- {u, (7 -|- H)~^v)). 
We will prove that the derivatives in r of both sides are equal, that is, 


-Trf(7 - 77)-i|:77) = Trf(7 + 77)-i|:77) + 




{u,fYi+H)-y 

{-lY + {u, {I+ H)-^v)' 


(3.12) 


(3.13) 


where we used the fact that ^ log (det (A)) = Tr (A A) if A is a square matrix depending 
smoothly on r. As a consequence, the two sides of (3.12) differ at most by a constant. But, 
since 77 —)• 0 as r —oo, both sides of (3.12) go to 0 as r —>■ oo, so the two sides are equal. 
Therefore our proof will be ready once we show that (3.13) holds. 

Since (7 — 77)“^ + (-^ + 77)“^ = 2(7 — 77^)“^, (3.13) is equivalent to 
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At this stage we use Lemma 3.2(iv) and then the cyclicity of the trace to obtain 
-2Tr(^(/ - = -2Tr(^(/ - = -2Tt:(^QH{I - 

Now note that if A is an x real symmetric matrix then Tt{QA) = —2r ^ij — 

—2r{u,Au), and thus 


Ar{u, = -2 Tv[{I-HY^^H). 


t 2\-1 d 


(3.15) 


On the other hand, on the right-hand side of (3.14) we may apply Lemma 3.2(v) and use 
the simple identity Qv = {—l)^Aru to get 


{u, |:(/ + H)-^v) = {u, (/ - HY HQv + 4r(/ - {Hu <S> u){I + H)-^v) 

= 4r{u, {I - H^)-^Hu) [(-1)^ + {u, {I + H)-\) 

Using (3.15) in (3.16) we get (3.14), which finishes the proof. 


(3.16) 

□ 


Proof of Lemma 3.2. 

(i) Fix A^ G N and r > 0, and define a upper triangular matrix S G as follows: 

c___ N-l-k ( 

~ ^j{ j-i with Ck — r (Af-i)! J 

for i,j G {0 ,... ,N — 1}. We claim that S is invertible, with inverse given by 

To check this, note first that, since both S and (as given above) are upper triangular, 
we have {SS~^)ij = 0 for f > j, while {SS~^)ii = SuSH = 1. Thus it remains to show 

that {SS-% = 0 when i < j. But {83-% = (-1)^^-^+^ Ei=i = 

I]fcZo(-l)^ k\{j-i-ky. ’ binomial theorem the last sum on the 

right-hand side is simply (—1 -|- 1)-^“* = 0. 

Now Lemma A.2 in Appendix A allows us to rewrite the symmetric matrix H in terms 
of Laguerre functions: for j > i, 

j 


Hi. = H., = 


k=i 


J - * 
k — i 


\{-l)%{2r^). 


(3.17) 


We will use this representation to show that S ^HS = H. We have 


{S-^HSfj = 


E 


k=i,...,N—l 


k — i 


(_l)2iV-2+i+fc^^^^ 


Note that the value of depends only on A; Letting ifn = V'n-i(2r^) — 'ifn{‘2.r‘^), so 
that Hki = (—l)'^V^fc-|-£_7v-i-i; and recalling that, by convention, = 0 for n < 0, we may 
write 

i ^ 


{s-^hs)„ = ^J2 


n=0 


E 


k=i,...,N—l, £=0,.--,i 
k-\-£—N-\-l=n 




k — i 


A^- 1 - 


1 ^^^_iy+k+N 


lAn. (3.18) 


Performing the change of variables k i—)• k + i, i i—>■ Al — 1 —and introducing the convention 
that {{{J = 0ifm>n>0, the sum in the square brackets turns into 

^N-l-P 


E 


k>0, 0<£<N-1 
k—l=n—i 


k 
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We claim now that the sum inO<^<A^ — 1 can be extended to ^ > 0. In fact, we may 
assume that k < N — 1 — i, since otherwise the hrst binomial coefficient vanishes. Since ^ 
is constrained to be ^ = A: + z — n<A^ — 1 — n<iV — 1, adding the terms with i > N does 
not really contribute to the sum. In view of this, and using Lemma A. I from Appendix A, 
our sum can be rewritten as 


E 


k>0, £>0 
£—k=i—n 


N-l-i 
k 


N -1- 




for i > n, 
for i < n. 


Now we substitute this formula into (3.18) and consider three separate cases: 

• If i = j, then {S-^HS)ii = Ylh=o (V) 

. If i < j, then (S-^HS)., = | | EL. (Ll) (-l)"V’n(2r2), 

where the second identity follows by summation by parts. 

• If i > j, then proceeding as for i < j we get {S~^HS)ij = ELo 

Applying Lemma A.3 from Appendix A we deduce that the last sum equals 

llw E {r-l)(-ir'/>«(2>'") = t E {r4)(-i)”'/'«(2r^). 

j<n<i 

In each case, the expression for {S~^HS)ij coincides with the formula for Hij in (3.17), 
which completes the proof of (i). 

(ii) We will use the contour integral representation of the Laguerre function 'tpn{x), 


ipnix) = e 


- ^-^/2 


1 


dt 


_ xt 

e i-‘ 


2m J t”+i(l-t)’ 


(3.19) 


where the integration is along a small circle around the origin (note that by Cauchy’s 
theorem this formula is consistent with our convention = 0 for n < 0). Together with 
the dehnition of L, (3.19) leads to 


POO 

J2r^ 


Ljk = I _ dx'iljj{x)'ijjkix) = j ^ j’j’ 


du dv 


XU _ XV 

1 — U 1 — 1) 


L// 


du dv 


ud+^il — u)u^+^(l — v) 

g \ 1 — u 1 — v J 


(3.20) 


(27ri)^ ff — uv) 

On the other hand, from the definition of H we get 


AT-l 


{H'^)jk = (-1)^^ ^ (?/)j+n_Ar(2r^) - l/)j+„_Ar-Hl(2r^)) (V’n-hfc-Ar(2r^) - i^n+k-N+ii'ir^)) 

n=0 

1 ff 

p® du dv 


(27ri)2 J J 


jV-l _2r2-2rhi,_ 2r^ 

g 1—u 1—v 


( 1 

1 ') 

( 1 

1 ') 

1 y^j + n—N+1 

y^j+n — N-\-2 J 

1 yU + k — +1 

yn-\-k — I\ -\-2 j 


1 


du dv 


n=0 

e 


-^'■^(l+i^ + i^)(l - (uu)^) 


(27ri)2 JJ — uv) 

The difference between (3.20) and (3.21) is then given by 


(3.21) 


Ajfc {H )jk — 


1 


(27ri 


D2 


du dv 


-2r^(l+^ + ^) , nAT 

o V ^ l-v J {^uv) 

+ — uu) 


Since 0 < j. A: < A^ — 1, the integrand has no poles in u and v inside the chosen contours, 
and hence the whole integral vanishes. 
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(iii) For the first formula, we compute directly Hu io get 


N-l 


{Hu)i = (-1)^^ ^ ^ ['0i+fc_Ar(2r2) _ = in{2r‘^) - = {Ri)i, 


k=0 


where the last identity follows because '0j_Ar(2r^) = 0 (since i < N). For the second 
one, we use the property -^{Ln(x) — Ln+i{x)) = Ln{x) of Laguerre polynomials to obtain 
^ (V’n(T) — 'ijjn+i{x)) = ^ (V’n(T) + il^n+iix)), which, together with the fact that -L„(0) = 1 
for all n € N, gives 

1 

- dx [i^nix) + ll^n+lix)] = '0n(2r^) - V'n-rl(2r^) 

^ Jo 


for all n G N. Hence we can write the entries of H as 

0 for f + j < — 1, 


Hij — < 


(-1)' 


for i + j = — 1, 


(3.22) 


{^i+j-N{‘Jr‘^) + '^i+j-N+i{‘Jr‘^)) for i + j > A^ - 1, 
with '&n(s) = jQdx'4>n{x) (note that ipo^x) = Now we can compute 


N-l 


N-l 


k=0 


k=N-i 

-(^'fc(2r2) + ^fc+i(2r2)) 


2 — 1 


MK-a-i)" 


= 2(-l)‘(l-e-' )-2^(-l) 

k=Q 

= 2(-l)'(l - e-^') - (-l)''I/o(2r2) + m,{2r^) = {R^y. 


(iv) From (3.22) we get 


Arr - 

dr^^j ~ 


lAf o—— 


This 


for i + j < A^ — 1, 
(—l)'^2re“^^ forz + j = A^ —1, 

y-l)^2r(ijji+j_N{‘^r‘^) + i’i+j-N+ii'^r'^)) ioi i + j > N - 1. 


\ \ / \ • I J \ / ' ' J IV // 

expression coincides with the i,j entry, for all i,j G {0,..., A^ — 1}, of the matrix 
QH^ which is given by 


2—1 

{QH)ij = -Ar'Y^Hkj - 2rHij = (-l)^2r('i/ij+j-Ar(2r^) + i/’i-ri-N-ri 
jb—n 


(2r2)) 


(v) For i, j G {0,..., A^ — 1} we have 

i-l N-l 

{QHyj + {HQ)ij = -4r(^Hij + Y,Hkj+ Y. 

k=0 k=j+l 


Since ^kj = J2k=o ^ik, the right-hand side of the last identity equals —4r J2k=o ^ik-, 

which coincides with —[ArHu ® u)ij. Thus, recalling the notation E = ArHu (8> u, we have 

QH = -HQ - E. (3.23) 


Now recall that if H is a square matrix which depends smoothly on a parameter r, then 
drA~^ = —A~^drAA~^. Then, in view of (iv) and the last identity, we have 

dr{I + H)-^ = -{I + H)-^QH{I + H)-^ = {I- hY\I - H){HQ + E){I + H)-^. 
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Comparing this with the right-hand side of the identity we seek to prove, we see that it 
is enough to check that HQ{I -|- H) + E = {I — H){HQ + E), which follows easily from 
(3.23). □ 


Appendix A. Some formulas for Hermite and Laguerre polynomials 

We begin with a combinatorial result which was used in the proof of Lemma 3.2(i) and 
which will also be used later in this appendix. 

Throughout this appendix we adopt the convention that, for A; G N and G Z, (^) = 0 if 
£ < 0 oi i > k (this can be justified, for instance, by replacing the factorials with Gamma 
functions). 


Lemma A.l. Let n, m G N and a G Z. Then 


E 

^,j>0 

j-i=a 


I J 
m 


(-ir = 


(-17 

(-1) 


m — n 


for a> 0, 


m-a 

. l{n>m-a} for O < 0. 

n — m + a J ^ ‘ 


Proof. Assume hrst that a > 0. Then the formula we seek to prove can be rewritten as 

a 


n\ z + a\ 

\ i ] \ m ] \m — n 

2>0 ^ ^ ^ ^ 


(A.l) 


For X G M and with our convention, we have (using Newton’s generalized binomial theorem) 


E E " 7“ = E ") +-) 

m>0 i>0 ^ ^ ^ '' i>0 ^ 


i-\-a 


= (1 +x)“ (-1 + (1 + x)r = (1 +x) V = 


By equating the coefficient in front of x™ on both sides, we obtain (A.l) . 
When o < 0, we first let 6 = —a > 0 and rewrite the desired identity as 


i>o 


+ ^ _h]hn>m+b}- 


n — m — 1 
n — m — h 


(A.2) 


Pick X G M such that |x| < 1 and |-r^| < 1. Using three times the identity 


(l_x)fc+i ^ \k 
^ ' n>k 


E 


n 


(A.3) 


(which is a straightforward consequence of Newton’s generalized binomial theorem) to¬ 
gether with our convention we have 


EE 

n>0 j>0 


n 


j + bj \mJ '' ' \jn 7(1 — x)l+^+^ 


(-l)l-™x^ = J]] 
i>o 

x^ {xjiX — x)) 




(-17- 


^b+m 


(l-x)^+i (l + x/(l-x))”^+^ (l-xY 


= X 


m+1 


E 

£> 6-1 


£ 

b-l 


X . 


(A.2) now follows from the fact that the coefficient of x” on the right-hand side is given 
by (YflYn-b) ''^1^6^ n>m + b and equals 0 when n < m + b. □ 
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Lemma A.2. For n, m E N with n > m and any r E M \ {0}, the following relation holds: 

'2™n!y/^ 

k=ra 


/ dx (Pnix)ipm{‘^r - x) = r™ ” 

Jr 

Similarly, for the case r = 0 we have 


2^m\ 


n — m 
k — m 


(-l)Vfe(2r2 


(A.4) 


I — ( 1) lm=T7. • 

Jr 


(A.5) 


Proof Recall that the Hermite polynomials have a simple generating function, namely 

Hn{x) _ ^2xt-P 


n=0 


n\ 


We write the convolution of Hermite functions in (A.4) as 

[ dx ipn{x)ipm(2r - x) = ^ == [ dxHn{x)e~''^/‘^Hm{‘ir - 

Jr v2”+"*7rn!m! Jr 

and then use the above generating function to evaluate the sum 


r7'7r>' /tr> 


n,m=0 


n!m! 


-2tit2 


Jr 

By equating the coefficient of tftf^ on each side, we obtain an explicit formula for the 
left-hand side of (A.4): 

/ dx Lpn{x)ipm{2r - x) 

Jr 


1 


\/2"'+”"n!m! 

1 


-P gn Q-m ^2r{ti+t2)-2tit2 


ti=t2=0 


y/2'^+mn\m\ 


£(-2)'f! 


£=0 


n\ m 


i V ^ 


(2r) 


n+m—2£ 


(A.6) 


In particular, we get (A.5), so from now on we will assume r ^ 0. 

Turning to the right-hand side of (A.4), we use the explicit power series expansion of 
the Laguerre polynomials, 


Lk{x) = 


e=o 


(-1)^ e 

’ X , 


£! 


(A.7) 


to rewrite it as 


2™n! 

2"m! 


1/2 fT. k 


EE 


n — m 
k — m 




-r2 {k\ {-'^Y (^2 


I £\ 


{2r^r- 


k=m £=0 

We need to show that this expression equals the right-hand side of (A.6) or, equivalently, 
that 


n—m k-\-m 

EE 

fc=o e=o 


n — m\fk + m\, , xfc_i_m. (—2r^)^ 


k 


I r ’ l\ 


E 


n-ir ’ il '■ 


where we have performed the changes of variables fe e->• fe -|- m on the left-hand side and 
i n — i on the right-hand side. Using our convention, this is equivalent to 


E 


£=0 lk=0 


E 


n — m\ fk + m 
k 


(- 1 ) 


n—ra—k 


{—2r'^Y ( m \ {—2r'^Y 


£\ 


e=o 


n — 


l\ 


It follows from Lemma A.l (or, more specihcally, from (A.l)) that the coefficients of 
on both sides of this identity coincide, and this finshes the proof. □ 
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Lemma A.3. For any n, m G N with n> m and any x G M, the following relation holds: 




n\ 


E 

fc =0 


n — k — 1 


n — m — 




k=m 


n — m 
k — m 




(A.8) 


Proof. We will use (A.7) in order to extract the coefficients of in the polynomials ap¬ 
pearing on both sides. The coefficient of x^ on the left-hand side of (A.8) is clearly 0 if 
i < n, while for n < i < n + m it is given by 

n\{i — n)\ — m — — nl(£ — n)l \i — m 

where we have used a variant of Vandermonde’s identity which can be obtained by equating 
the coefficient of x” in the expansion of both sides of the identity = 

obtained by using (A.3). On the other hand, for ^ < m the coefficient of x^ on 
the right-hand side of (A.8) is clearly zero, while for m<i<n + m it is given by 



(^_lY m _ (-1)^+^ / m \ 

ml{i — m)\ \ k J \£ — m J ml(i — m)\ \£ — nj ’ 


(A.IO) 


where we used the change of variables k ^ k + m and Lemma A.l. Notice that (A.IO) 
equals 0 by our convention if m < .^ < n and it clearly equals (A.9) ifn<£<n-|-m (recall 
that we are assuming n > m). The proof is thus complete. □ 


Appendix B. Proof of Lemma 2.3 


Throughout the proof we will use ci and C2 to denote positive constants which do not 
depend on L and whose value may change from line to line. We will denote by || • ||i and 
II • II2 the trace class and Hilbert-Schmidt norms of operators on L^(M). We recall that 

P^lli < II-4II2IIBII2 and \\A\\l = j dxdyA{x,y)‘^ (B.l) 

if A has integral kernel A{x,y); for more details see [QR14, Section 2] or [Sim05]. 

In view of (2.11) we write where 

® ^lieTm,N‘rcosh{L)\^ ^[—L,0])^^V® ^[0,L]) ^^cosh(L)® ''Herm,Af) 

^ ''Herm,Af(C L,0]W^V® ^[ 0 ,L] W^cosh(L)^ ''Herm,Af- 


(d1) 


We will focus on Q 
very similar so we will omit it. 

We factor as 


^ and show that it goes to zero in trace norm, the proof for 


(d2) . 


IS 




with Ti = 


'^Herm,A^r 7-cos' 
(r) 


;h(L) (' 


,-LD 


_ p(^) 'ip 

LDi 


-L,0]> 


and T2 = Pr(e-^ - R[o,i]) Prcosh(L)e KHerm,Af- 


By (B.l), it is enough to show that IIT1II2IIT2II2 —^ 0 as L —)• 00. We start with Ti, 
which is made of two terms which we will bound separately. By (B.l) and the fact that 
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the family ((/?n)neN is orthonormal we have 

poo pr 

||e'^“KHerm,NPrcosh(L)e"^°Pr||i = i dx dy ^Pnix)(Pn'(x) 

„l-r,J-00 J-QO 


n,n'=0 ' 
X 


i 


rcosh(L),oo)^ 


dzdz' (Pn{z)(fn'{z^)e ^^iz,y)e 


TVl / /* 

■i!' / 

n=0 


dz(pn{z)e ^^{z,y) 


rcosh(L) 

/*oo /*oo 

«/‘rcosh(L) J —oo 

where we have used the Cauchy-Schwarz inequality. Using the formula for the kernel of 
which is implicit in (2.6) and (2.7) we see that that the y integral is just a Gaussian 
integral, and computing it gives 

roo 

||e^^KHermjvPrcosh(L)e~^^Pr||i < coth(L)-l / ^2L-zH^nh(L) _ 

II Herm,A rcosh(L) r||2- Jrcosh{L) 

The last integral is bounded by for all L > 0, and thus, recalling 

that we are assuming r > 0, 

W LDw p —LDp ||2 ^ _ „2NL—C2e^^ 

||6 ■'HemijAfr rcosh(L)® i^r’ll2 — *^1® 


for sufficiently large L. The estimate for the other term appearing in Ti is very similar 
and leads to the same type of bound. We deduce that 


IT1II2 < cie 


NL-C2e^^ 


(B.2) 


for large enough L. On the other hand it is easy to check that the same calculation as 
above leads to 

IIT2II2 <cie^^ (B.3) 

(note that in this case the projection P^coshjL) appearing in Ti is replaced by Pr-cosh(L); 
this accounts for the fact that the factor disappears from the upper bound). By 

combining (B.2) and (B.3) together we immediately get 

->0. 

L—¥ 00 
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